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Abstract. We derive the astrometric orbit of the photo- 
center of the close pair a UMi AP (= a UMi Aa) of the 
Polaris multiple stellar system. The orbit is based on the 
spectroscopic orbit of the Cepheid a UMi A (orbital pe- 
riod of AP: 29.59 years), and on the difference A/i between 
the quasi-instantaneously measured HIPPARCOS proper 
motion of Polaris and the long-term-averaged proper mo- 
tion given by the FK5. There remains an ambiguity in 
the inclination i of the orbit, since A/x cannot distinguish 
between a prograde orbit (i = 50 ?1) and a retrograde 
one (i = 130 ?2). Available photographic observations of 
Polaris favour strongly the retrograde orbit. For the semi- 
major axis of the photo-center of AP we find about 29 
milliarcsec (mas). For the component P, we estimate a 
mass of 1.5 Mq and a magnitude difference with respect 
to the Cepheid of 6.5 mag. The present separation between 
A and P should be about 160 mas. 

We obtain the proper motion of the center-of-mass of 
a UMi AP with a mean error of about 0.45 mas/year. Us- 
ing the derived astrometric orbit, we find the position of 
the center-of-mass at the epoch 1991.31 with an accuracy 
of about 3.0 mas. Our ephemerides for the orbital correc- 
tion, required for going from the position of the center- 
of-mass to the instantaneous position of the photo-center 
of AP at an arbitrary epoch, have a typical uncertainty 
of 5 mas. For epochs which differ from the HIPPARCOS 
epoch by more than a few years, a prediction for the actual 
position of Polaris based on our results should be signif- 
icantly more accurate than using the HIPPARCOS data 
in a linear prediction, since the HIPPARCOS proper mo- 
tion contains the instantaneous orbital motion of about 
4.9 mas/year = 3.1 km/s. Finally we derive the galactic 
space motion of Polaris. 

Key words: astrometry - binaries: general - Cephcids 
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Fig. 1. Astrometric orbit (prograde or retrograde) of the 
photo-center of a UMi AP. The retrograde orbit is our pre- 
ferred solution. For detailed explanations see Sect. 3.2.4. 

1. Introduction 

Polaris (a Ursae Minoris, HR 424, HD 8890, ADS 1477, 
FK 907, HIP 11767) is a very interesting and important 
object, both from the astrophysical point of view and from 
the astrometric one. For astrophysics, the most remarkable 
feature of the multiple stellar system Polaris is the fact 
that its main component, namely a UMi A, is a Cepheid 
variable with a very unusual behaviour. In astrometry, Po- 
laris is one of the most frequently and accurately observed 
objects, mainly because it is located so close to the North 
celestial pole and can be used for calibration purposes. 

Up to now, the binary nature of Polaris was essentially 
neglected in ground-based fundamental astrometry, e.g in 
the FK5 (Fricke et al. 1988). This was justified by the 
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limited accuracy reached by the meridian-circle observa- 
tions. Now, the high-precision astrometric measurements 
carried out with the HIPPARCOS satellite (ESA 1997) 
require strongly to take into account the binary nature 
of Polaris in order to obtain an adequate astrometric de- 
scription of a UMi. Similar procedures are required for 
many other binaries among the fundamental stars in or- 
der to be included properly into the Sixth Catalogue of 
Fundamental Stars (FK6; Part I: Wielen ct al. 1999c; see 
also Wielen et al. 1998). 

The main purpose of the present paper is to obtain 
a reliable astrometric orbit for Polaris, and to use this 
astrometric orbit for obtaining high-precision values for 
the position and proper motion of Polaris. This is done by 
combining the known spectroscopic orbit of Polaris with 
ground-based astrometric data given in the FK5 and with 
the HIPPARCOS results. Before doing so in the Sects. 3 
and 4, we present in Sect. 2 an overview of the Polaris 
system. 

2. Overview of the Polaris System 

Polaris is a multiple stellar system, which consists of a 
close pair, a UMi A and a UMi P (= a UMi a), and a 
distant companion, a UMi B, and two distant components, 
a UMi C and a UMi D. We use here the designation 'P' 
for a close companion of A, which was used in the IDS 
and was adopted by the CCDM and by the HIPPARCOS 
Input Catalogue, rather than the traditional version 'a', 
which is used e.g. by the WDS and by CHARA). 

2.1. The Cepheid a UMi A 

The main component of Polaris is a low-amplitude 
Cepheid with a pulsational period of about 3.97 days. 
This period is increasing with time (e.g. Kamper & Fer- 
nic 1998). According to Feast & Catchpole (1997), a UMi 
A is a first-overtone pulsator (rather than a fundamen- 
tal one), since a UMi A is too luminous for a fundamen- 
tal pulsator, if they apply their period-luminosity relation 
(for fundamental pulsators) to Polaris. The fundamental 
period of a UMi A would follow as Pq = 5.64 days, if 
the observed period is the first-overtone period P\ (using 
the relation between P\ and Po derived by Alcock et al. 
(1995) for Galactic Cephcids). An extraordinary property 
of a UMi A among the Cepheids is that the amplitude 
of its pulsation has been dramatically declined during the 
past 100 years, as seen both in the light curve and in 
the radial-velocity curve (Arellano Ferro 1983, Kamper & 
Fcrnic 1998, and other references given therein). The full 
amplitude was about 07*12 in my and about 6 km/s in 
radial velocity before 1900, and seems now to be rather 
constant at a level of only OTTO in my and at 1.6 km/s in 
radial velocity. An earlier prediction (Fernie et al. 1993) 
that the pulsation should cease totally in the 1990s was 
invalid. A discussion of the HIPPARCOS parallax and of 



the absolute magnitude of a UMi A is given in the next 
Sect. 2.2. 

2.2. The spectroscopic- astrometric binary a UMi AP 

The Cepheid a UMi A is a member of the close binary 
system a UMi AP. This duplicity was first found from the 
corresponding variations in the radial velocity of a UMi 
A. However, the interpretation of the radial velocities of 
a UMi A in terms of a spectroscopic binary is obviously 
complicated by the fact that a UMi A itself is pulsating 
and that this pulsation varies with time. We use in this 
paper the spectroscopic orbit derived by Kamper (1996), 
which is based on radial velocity observations from 1896 
to 1995. Kamper (1996) took into account changes in the 
amplitude of the pulsation and in the period of pulsation, 
but used otherwise a fixed sinusoid for fitting the pulsa- 
tion curve. In an earlier paper, Roemer (1965) considered 
even 'annual' changes in the form of the pulsation curve. 
In Table 4, we list the elements of the spectroscopic orbit 
of A in the pair AP given by Kamper (1996, his Table III, 
DDO + Lick Data). The orbital period of a UMi AP is 
29.59 ±0.02 years, and the semi-amplitude is K a = 3.72 
km/s. The value of aAsinj = 2.934 AU corresponds to 
about 22 milliarcsec (mas), using the HIPPARCOS paral- 
lax. 

Attempts to observe the secondary component a UMi 
P directly or in the integrated spectrum of a UMi AP 
have failed up to now. Burnham (1894) examined Polaris 
in 1889 with the 36-inch Lick refractor and found no close 
companion to a UMi A (nor to a UMi B). Wilson (1937) 
claimed to have observed a close companion by means of 
an interferometer attached to the 18-inch refractor of the 
Flower Observatory Jcffcrs (according to Roemer (1965) 
and to the WDS Catalogue) was unable to confirm such a 
companion with an interferometer at the 36-inch refractor 
of the Lick Observatory. HIPPARCOS (ESA 1997) has not 
given any indication for the duplicity of Polaris. Speckle 
observations were also unsuccessful (McAlister 1978). All 
these failures to detect a UMi P directly are not astonish- 
ing in view of the probable magnitude difference of A and 
P of more than 6 m and a separation of A and P of less 
than 0'.'2 (see Sect. 3.2.5). Roemer and Herbig (Roemer 
1965) and Evans (1988) searched without success for light 
from a UMi P in the combined spectrum of a UMi AP. 
From IUE spectra, Evans (1988) concluded that a main- 
sequence companion must be later than A8V. This is in 
agreement with our results for a UMi P, given in Table 5. 
A white-dwarf companion is ruled out by the upper limit 
on its effective temperature derived from IUE spectra and 
by considerations on its cooling age, which would be much 
higher than the age of the Cepheid a UMi A (Landsman 
et al. 1996). 

After Polaris had become known as a long-period spec- 
troscopic binary (Moore 1929), various attempts have 
been made to obtain an astrometric orbit for the pair a 
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UMi AP. Meridian-circle observations were discussed by 
Gerasimovic (1936) and van Herk (1939). While van Hcrk 
did not find a regular variation with a period of 30 years, 
Gerasimovic claimed to have found such a modulation. 
However, the astrometric orbit of the visual photo-center 
of a UMi AP determined by Gerasimovich (1936) is most 
probably spurious, since he found for the semi-major axis 
of the orbit a p h(AP) ~ 110 mas, which is much too high 
in view of our present knowledge (a ph ( AP ) = 29 mas). 
More recent meridian-circle observations gave no indica- 
tions of any significant perturbation. This is not aston- 
ishing in view of the small orbital displacements of the 
photo-center of AP of always less than 0'.'04. Long-focus 
photographic observations have been carried out at the Al- 
legheny Observatory (during 1922-1964), the Greenwich 
Observatory, and the Sproul Observatory (during 1926- 
1956), mainly with the aim to determine the parallax of 
Polaris. The discussion of this material by Wyller (1957, 
Sproul data) and by Roemer (1965, Allegheny data) did 
not produce any significant results. The Allegheny plates 
were later remeasured and rediscussed by Kamper (1996), 
using his new spectroscopic orbital elements. Kamper also 
rediscussed the Sproul plates. While the Sproul data gave 
no relevant results for a UMi AP, the Allegheny data gave 
just barely significant results, such as a p h(AP) = 19.5 ±6.5 
mas. For our purpose (see Sect. 3.2.3), the most important 
implication derived by Kamper (1996) from the Allegheny 
data is that the astrometric orbit of AP is most probably 
retrograde, not prograde. 

In Sect. 3 we shall present a more reliable astromet- 
ric orbit of a UMi AP by combining ground-based FK5 
data with HIPPARCOS results, using Kamper's (1996) 
spectroscopic orbit as a basis. 

The HIPPARCOS astrometric satellite has obtained 
for a UMi AP a trigonometric parallax of pn = 7.56 ±0.48 
mas, which corresponds to a distance from the Sun of 
r H = 132 ± 8 pc. In the data reduction for HIPPARCOS, 
it was implicitely assumed that the photo-center of the 
pair AP moves linearly in space and time, i.e. a 'standard 
solution' was adopted. This is a fairly valid assumption, 
since the deviations from a linear fit over the period of ob- 
servations by HIPPARCOS, about 3 years, are less than 1 
mas (see Sect. 4.2). Hence the HIPPARCOS parallax ob- 
tained is most probably not significantly affected by the 
curvature of the orbit of AP. Nevertheless, it may be reas- 
suring to repeat the data reduction of HIPPARCOS for a 
UMi, adopting the astrometric orbit derived here for im- 
plementing the curvature of the orbit of the photo-center 
of a UMi AP. 

The mean apparent visual magnitude of the combined 
components A and P is my.AP = 1-982 (Feast & Catch- 
pole 1997). This agrees fairly well with the HIPPARCOS 
result (ESA 1997) my,AP = 1-97. In accordance with most 
authors we assume that the reddening Eb-v and the ex- 
tinction Ay of the Polaris system are essentially zero ( e.g., 
Turner 1977, Gauthicr & Fcrnie 1978), within a margin of 



±0.02 in E B -y and ±0.06 in Ay. Using the HIPPAR- 
COS parallax, we find for the mean absolute magnitude 
of AP M v ,ap = -3.63 ± 0.14. If we use our results of 
Table 5 for component P, i.e. My,p ~ ±2.9, and subtract 
the light of P from Mv,ap, then the absolute magnitude 
of the Cepheid component A is My, a = —3.62 ± 0.14. 
Unfortunately, the pecularities in the pulsation of a UMi 
A are certainly not very favourable for using this nearest 
Cepheid as the main calibrator of the zero-point of the 
period-luminosity relation of classical Cepheids. 

2.3. The visual binary a UMi (AP)-B 

Already in 1779, W. Herschel (1782) discovered the visual- 
binary nature of Polaris. The present separation between 
AP and B is about 18'.' 2. This separation corresponds to 
2400 AU or 0.012 pc, if B has the same parallax as AP. 
Kamper (1996) has determined the tangential and radial 
velocity of B relative to AP. Both velocities of B agree 
with those of AP within about 1 km/s. Hence Kamper 
(1996) concludes that B is most probably a physical com- 
panion of AP, and not an optical component. The phys- 
ical association between AP and B is also supported by 
the fair agreement between the HIPPARCOS parallax of 
AP (rn = 132 ± 8 pc) and the spectroscopic parallax of B 
(114 pc, as mentioned below). 

The spectral type of B is F3V. The magnitude differ- 
ence between B and the combined light of AP is Amy = 
6.61 ±0.04 (Kamper 1996). Using m v ,AP = 1-98, this im- 
plies for B an apparent magnitude of my,B = 8.59 ± 0.04. 
Adopting the HIPPARCOS parallax (and no extinction), 
we obtain for B an absolute magnitude of Mv,b = ±2.98± 
0.15. The standard value of My for an F3V star on the 
zero-age main sequence is ±3.3. If we use this standard 
value for My, we obtain for B a spectroscopic distance 
of r = 114 pc. Similar values of the spectroscopic dis- 
tance were derived (or implied) by Fernie (1966), Turner 
(1977), and Gauthier & Fernie (1978). These authors were 
interested in the absolute magnitude (and hence in the 
distance) of B in order to calibrate the absolute magni- 
tude of the Cepheid A. Now the use of the HIPPARCOS 
trigonometric parallax is, of course, better suited for this 
purpose. 

The typical mass of an F3V star is Mb = 1.5 Mq. If 
we use for the masses of A and P the values adopted in 
Table 5 (6.0 ± 1.54 Mq), we obtain for the triple system a 
total mass of Mtot = 9.0 Mq). We derive from Pb-ap = 
18'.' 2 and the statistical relation a = 1.13 p an estimate 
for the semi-major axis of the orbit of B relative to AP of 
ob-ap ~ 21" or 2700 AU. From Kepler's Third Law, we 
get then an estimate of the orbital period of B, namely 
Pb ~ 50 000 years. 

From the data given above, we can estimate the ac- 
celeration <7ap °f the center-of-mass of the pair a UMi 
AP due to the gravitational attraction of a UMi B. If we 
project this estimate of c/ap on one arbitrarly chosen di- 
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rection, we get for AP a typical 'one-dimensional' acceler- 
ation of about 0.003 (km/s) /century or 0.4 mas/century 2 . 
Therefore, we should expect neither in the radial velocity 
nor in the tangential motion of AP a significant deviation 
from linear motion due to the gravitational force of B dur- 
ing the relevant periods of the observations used. For all 
present purposes, it is fully adequate to assume that the 
ccnter-of-mass of the pair a UMi AP moves linearly in 
space and time. The same is true for the motion of B. 

A modulation of the relative position of B with respect 
to the photo-center of AP with a period of about 30 years 
is not seen in the available observations of B. This is in 
accordance with our determination of the motion of the 
photo-center of AP with respect to the cms of AP, given in 
Table 7. The expected amplitude of the modulation is less 
than 0'.'04 and is obviously not large enough with respect 
to the typical measuring errors in the relative position of 
B. 

The contribution of the orbital motion of the center-of- 
mass (cms) of AP, due to B, to the total space velocity of 
AP is of the order of a few tenth of a km/s. The expected 
value of the velocity of B relative to the cms of AP is of 
the order of 1 km/s. 

2.4. a UMi C and a UMi D 

In 1884 and 1890, Burnham (1894) measured two faint 
stars in the neighbourhood of a UMi AB. In 1890.79, the 
component C had a separation of 44 '.'68 from A, and the 
component D 82 '.'83. According to the WDS Catalogue, 
the apparent magnitudes of C and D are 13!™1 and 127*1. 

The nature of the components C and D is unclear. 
The probability to find by chance a field star of the corre- 
sponding magnitude with the observed separation around 
a UMi A (galactic latitude b = +26?46 (Wielen 1974)) 
is of the order of 10 percent for each component. This 
favours on statistical grounds a physical relationship of 
the components C and D with A. If C and D are physical 
members of the Polaris system (instead of being optical 
components), their absolute magnitudes in V would be 
+ 7" l 5 and + 6!™5. Due to the low age of the Polaris sys- 
tem of about 70 million years (deduced from the Cephcid 
a UMi A), they would either just have reached the zero- 
age main sequence, or they may still be slightly above this 
sequence (i.e. pre-main-sequence objects, Fernie 1966). 

3. Astrometric orbit of a UMi AP 

In this section we determine the astrometric orbit of the 
photo-center of the pair a UMi AP (i.e. essentially of A) 
with respect to the center-of-mass of AP. We adopt all 
the elements of the spectroscopic orbit of A in the system 
AP, derived by Kampcr (1996). The remaining elements, 
i.e. the orbital inclination i and the nodal length f2, are 
basically obtained from the following considerations: 



The observed difference A/x between the instantaneous 
proper motion of a UMi A, provided by HIPPARCOS for 
an epoch T Ci h ~ 1991.31, and the mean proper motion of 
a UMi A, provided by long-term, ground-based observa- 
tions, summarized in the FK5, is equal to the tangential 
component of the orbital velocity of A with respect to the 
center-of-mass of the pair AP. Using the spectroscopic or- 
bit of A and the HIPPARCOS parallax, we can predict 
A/i for various adopted values of i and ft. Comparing the 
predicted values of A/i with the observed difference A/i, 
we find i and f2. The length of the two-dimensional vector 
of A/i gives us the inclination i; the direction of Art fixes 
then the ascending node fi. Unfortunately, two values of 
i, namely i and 180° — i, predict the same value for A/i 
(see Fig. 1). This ambiguity corresponds to the fact that 
A/i itself does not allow us to differentiate between a pro- 
grade orbit and a retrograde one. In the case of a UMi 
AP, it is fortunate that the ground-based observations of 
the Allegheny Observatory strongly favour the retrograde 
orbit over the prograde one. 

3.1. The determination o/A/i 

3.1.1. The proper motion of the center-of-mass of AP 

We determine first the proper motion /i cms (AP) of the 
center-of-mass of the pair AP. (/i is used here for /x a * = 
/t Q cos<5 or for fj,g). The proper motion /ifks of a UMi 
given in the FK5 should be very close to /i cms (AP)j since 
the ground-based data are averaged in the FK5 over about 
two centuries, which is much larger than the orbital pe- 
riod of AP of about 30 years. In Table 1, we list /*fk5 
in the FK5 system and, by applying appropriate system- 
atic corrections, in the HIPPARCOS/ICRS system. The 
mean errors of /i F K5 in the HIPPARCOS system include 
both the random error of /xfks and the uncertainty of the 
systematic corrections. 

Another determination of /i cm s(AP) is based on the po- 
sitions £h(2c,h) and xfks(7c j fk5) at the central epochs 
T c ,h and T CiFK5 of the HIPPARCOS Catalogue and of 
the FK5. The designation x stands for a* = acos<5 or 
5, where a is the right ascension and 5 the declination 
of a UMi. The position £ F K5(Tc,FK5) represents a time- 
averaged, 'mean' position in the sense of Wielen (1997). 
Before being used, xfksCTc.fks) must be reduced to the 
HIPPARCOS/ICRS system. 

The HIPPARCOS position is (approximately) an 'in- 
stantaneously' measured position of the photo-center of 
AP. Before combining the HIPPARCOS position with 
xfk5 to a mean proper motion /io, we have to re- 
duce xh to the mean position a; meanp h(AP),H(T'c,H) of 
the photo-center of AP at time T c .h- This is done 
by going first from Xk(T c> k) to the center-of-mass 
x cms(AP).H(Tc,H) by subtracting from xh the orbital dis- 
placement Ax or b, p h(AP)(Tc,H) predicted by the astromet- 
ric orbit of the photo-center of AP. Then we have to add to 
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Prograde orbit Retrograde orbit 

(Preferred solution) 



Quantity [mas/year] 


System 




m.e. 


M<S 


m.e. 




m.e. 


Mi 


m.e. 


MFK5 


FK5 


+ 38.30 


0.23 


- 15.20 


0.35 


+ 38.30 


0.23 


- 15.20 


0.35 


systematic correction 




+ 3.20 


0.94 


- 1.53 


0.66 


+ 3.20 


0.94 


- 1.53 


0.66 


A*FK5 


HIP 


+ 41.50 


0.97 


- 16.73 


0.75 


+ 41.50 


0.97 


- 16.73 


0.75 


Mo 


HIP 


+ 41.05 


0.58 


- 15.05 


0.45 


+ 40.56 


0.58 


- 14.67 


0.45 


A*FK5 — MO 


HIP 


+ 0.45 


1.13 


- 1.68 


0.87 


+ 0.94 


1.13 


- 2.06 


0.87 


Mm 


HIP 


+ 41.17 


0.50 


- 15.49 


0.39 


+ 40.81 


0.50 


- 15.22 


0.39 


MH 


HIP 


+ 44.22 


0.47 


- 11.74 


0.55 


+ 44.22 


0.47 


- 11.74 


0.55 


Afl = mh - Mm 


HIP 


+ 3.05 


0.69 


+ 3.75 


0.67 


+ 3.41 


0.69 


+ 3.48 


0.67 



Table 2. Proper-motion difference Am between mh and Mm of the photo-center of a UMi AP at the epoch T c ,h = 
1991.31 









Prograde orbit 


Retrograde orbit 












(Preferred solution) 


Quantity 


Unit 


Observed 


Predicted 


Observed 


Predicted 


Afi at 


mas/year 


+ 3.05 


± 0.69 


+ 3.05 


+ 3.41 ± 0.69 


+ 3.41 


Ams 


mas/year 


+ 3.75 


± 0.67 


+ 3.75 


+ 3.48 ± 0.67 


+ 3.48 


A/xtot 


mas/year 


4.83 


± 0.61 


4.83 ± 0.32 


4.87 ± 0.61 


4.88 ± 0.32 


6a m 


o 


39.1 


± 8.8 


39.1 ± 3.5 


44.4 ± 8.8 


44.4 ± 3.4 



x cms(AP),H(^c,H) the (constant) off-set between the mean 
position of the photo-center x mC anph(AP) an d the center- 
of-mass (see Fig. 1). Using now a* and 5, we obtain 

Q;*,meanph(AP),H(7 1 Cl H) = a*,cms(AP),H(^c,H) 
3 

— - ea ph (AP)(cosa>sm£! + sinwcosficosi) , (1) 

^mean ph(AP),H(^ 1 c,H) — ^cms(AP),H(^ 1 c,H) 
3 

— - e a ph ( AP ) (cos uj cos — sin u sin £1 cos i) , (2) 

where a p h(AP) is the semi-major axis of the orbit of the 
photo-center of AP around the center-of-mass of AP. The 
other elements of this orbit are: eccentricity e, inclination 
i, longitude of periastron lj, position angle of the ascend- 
ing node fi, orbital period P, epoch of periastron passage 
T per i. The quantities in the Eqs. (1) and (2) which follow 
after — | e are just the Thiele-Innes elements B and A. 
The equations use the fact that, in the orbital plane, the 
time-averaged position is located on the major axis, to- 
wards the apastron, at a distance of § e a from the center- 
of-mass. 

If a* and 8 would change linearly with time, we could 
determine the mean proper motion mo from 

2-mean ph(AP) ,H 



However, for Polaris we should use more accurate formulae 
because it is so close to the celestial pole. We determine 
Mo strictly by requiring that Mo(Pc,h) is that proper mo- 
tion which brings the object from £ me anph(AP),H(Tc,H) to 
£fk5(T Cj fk5)- For calculating the (small) foreshortening 
effect, we have adopted the radial velocity of the center- 
of-mass of AP, v T = 7 = —16.42 km/s (Kamper 1996). 

The agreement between the two mean proper motions 
Mfks and Mo is rather good (Table 1). For determining 
the best value Mm of the mean motion of the photo-center, 
which is equal to the proper motion of the center-of-mass, 
we take the weighted average of MFK5 and mo- Since the 
orbital corrections to ih are different for the prograde and 
retrograde orbits, we have two values for mo and hence for 
Mm- In both cases, we had to iterate the determinations of 
the orbital elements (i and fi) and of Mo (and hence /j m ), 
since Mo depends on the orbital corrections. The values for 
Mm = Mcms finally adopted are listed in Table 1. 

3.1.2. The HIPPARCOS proper motion mh 

The HIPPARCOS proper motion mh of Polaris (ESA 
1997) refers to the photo-center of AP. Basically, mh is 
the sum of the proper motion Mcms(AP) of the center-of- 
mass (cms) of AP and of the orbital motion AM or b,ph(AP) 
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(abbreviated as A/i) of the photo-center of AP with re- 
spect to the cms of AP at time T Cj h: 



MhCT C) h) = Mcms(AP) + A/i(T CiH ) • 



(4) 



During the reduction of the HIPPARCOS data, a linear 
'standard' solution was applied to Polaris. The variation 
of A/i(i) during the period of observations of about three 
years was neglected. This slightly complicates the com- 
parison of the observed A/i with the orbital ephemerides. 
In Sect. 3.1.5 we assume that /in is obtained from a 
linear fit to quasi-continuously measured true positions 
over a time interval £>h, centered at time T Cj h. From the 
correlation coefficients given in the HIPPARCOS Cata- 
logue, we derive for the central epochs T a ,n — 1991.26 
and T^h = 1991.35. We neglect the slight difference be- 
tween Tq, ; h and Tgn and use the average of both, namely 
T c , h = 1991.31. From the epochs of the individual observa- 
tions of Polaris by HIPPARCOS, we estimate D H = 3.10 
years. 



3.1.3. The observed value of A/i 

The observed value of A/i is derived from 

A/j(T CjH ) = Mh(Tc,h) - Mcms(AP)(r Ci H) ■ (5) 

The values of A/i in a* and S, derived from Eq. (5), 
are listed in Tables 1 and 2. Table 2 gives also the total 
length A/itot of the A/i vector, 



(A/xtot) 2 = (A/w) 2 + (A^) 2 



(6) 



and the position angle 6a^ of the A/i vector, 

Oa m = arctan(A/i Qst /A/i,5) . (7) 

All the values are valid for the equinox J2000 in the HIP- 
PARCOS/ICRS system at the epoch T C , H = 1991.31. 

Since /i cms depends on the direction of motion in 
the orbit (prograde or retrograde), this is also true for 
A/i, and we obtain therefore two values for A/i. Table 2 
shows that A/i to t ~ 5 mas/year and ©a^, are statistically 
quite significant and rather well determined. A value of 
A/itot = 4.87 mas/year corresponds to a tangential veloc- 
ity of 3.05 km/s. Hence the 'instantaneous' HIPPARCOS 
proper motion of Polaris has a significant 'cosmic error' 
(Wielen 1995a, b, 1997, Wielen et al. 1997, 1998, 1999a, b) 
with respect to the motion of the center-of-mass. If Polaris 
were not already known as a close binary, our A/i method 
(Wielen et al. 1999a) would have detected Polaris to be a 
A/i binary because of its large test parameter Ffh = 6.18 
for /i FK 5 - Mh- 

3.1.4. The photo-center of a UMi AP 

The HIPPARCOS observations refer to the photo-center 
of a UMi AP, since the pair is not resolved by HIP- 
PARCOS. The 'phase' used in constructing the HIPPAR- 
COS Catalogue is practically identical to the phase of the 



photo-center, because the magnitude difference Atoap of 
more than 6 mag between A and P is quite large and be- 
cause the separation between A and P at Th was rather 
moderate (about 93 mas). It can also be shown that the 
component B does not significantly affect the HIPPAR- 
COS measurements of AP, because of Atov,b-ap = 6.61, 
in spite of its separation p = 18". The HIPPARCOS ob- 
servations have been carried out in a broad photometric 
band called Hp. The photo-center refers therefore to this 
photometric system. 

The spectroscopic orbit, however, refers to component 
A. We have therefore to transform the value aAsini of 
the spectroscopic orbit into a p h(AP) sini for obtaining an 
astrometric orbit of the photo-center of AP. The relation 
between oa and a p h(AP) is given by 



ftph(AP) = (1 - -g) a A 



(8) 



where B and [3 are the fractions of the mass M and the 
luminosity L of the secondary component P: 



B 



P = 



Mi 



Ma + M p ' 

L P 1 



L A + Lp 1 + 100.4 Am AP ' 

AniAp is the magnitude difference between A and P: 



(9) 
(10) 



AmAP = rnp — toa = Mp — Ma ■ (11) 
Using the results given in Table 5, we find for a UMi AP 



1 - ^ = 0.988 , 
B 



(12) 



with an estimated error of about ±0.010. For calculating 
(3, we have assumed that AmAP is the same in Hp as in V. 
This approximation is fully justified for our purpose. We 
derive (see the end of Sect. 3.1.6) from Kampcr (1996) for 
component A: 

a A sin i = 2.934 ±0.028 AU . (13) 
The HIPPARCOS parallax of Polaris (ESA 1997) is 
p H = 7.56 ±0.48 mas . (14) 
This leads to 



aA sinz = 22.18 ± 1.42 mas . 



(15) 



Using Eqs. (8), (12), and (15), we obtain for the photo- 
center 



a ph (Ap) sinz = 21.91 ± 1.42 mas 



(16) 
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3.1.5. The predicted value of A/i 

For predicting A/i, we use four elements (P, e, T per i, a;) of 
the spectroscopic orbit derived by Kamper (1996), and 
a ph(AP) sm « according to Eq. (16), all listed in Table 4 . 
In addition we adopt various values of i and Q in order to 
produce predicted values of A/i at time T C) h as a function 
of i and ft. 

Since the observed value of A/i is not an instanta- 
neously measured tangential velocity, we mimic the HIP- 
PARCOS procedure of determining /in- We calculate the 
positions Ax or bit jP h(AP)(^) = Ax(i) of the photo-center of 
AP with respect to the cms of AP as a function of time, 
using standard programs for the ephemerides of double 
stars. We then carry out a linear least-square fit to these 
positions over a time interval of length Djj = 3.10 years, 
centered at T c H = 1991.31: 



+£>h/2 



Ax av (T c . H ) = / Ax(T C! h + r)dT 

L>H J 



--Dh/2 
+£>h/2 



A/i av (r C:H ) 



12 f 
M J 



A:e(T C! h + t)t dr . 



(17) 



(18) 



Table 3. Determination of the inclination 


i from A/itot 


Prograde orbit 


Retrograde orbit 






(Preferred solution) 




A/itot 




A/itot 




[mas /years] 




[mas /years] 


Observed: 


4.83 ± 0.61 


Observed: 


4.87 ± 0.61 


Predicted 




Predicted 




for i [°]: 




for i ["]: 




30 


9.32 


150 


9.32 


40 


6.58 


140 


6.58 


45 


5.63 


135 


5.63 


50 


4.85 


130 


4.85 


55 


4.19 


125 


4.19 


60 


3.64 


120 


3.64 


70 


2.78 


110 


2.78 


50.1 


4.83 ± 0.32 


130.2 


4.88 ± 0.32 



-D H /2 



Tests have shown that especially A/i av is not very sensitive 
against small changes in the slightly uncertain quantity 
Dr. Actual numbers for the predicted values A/i av (T C; H) 
are given in the Tables 2 and 3. 

3.1.6. The problem o/T per i and of a a sin i 

In his paper, Kamper (1996, his Table III) gives for his 
best orbit (DDO+Lick Data) a value for T pcri = 1928.48 ± 
0.08. This is exactly the value derived by Roemer (1965) 
from the Lick Data and also quoted in Kamper's Table 
III under 'Lick Data'. There are three possibilities for this 



For our purpose, a value of T per i closer to T Ci H should 
be chosen. Using T por i = 1928.48 ± 0.12 and P = 29.59 ± 
0.02 years, we obtain an alternative value (two periods 
later) of 



Tp eri = 1987.66 ±0.13 . 



(19) 



We have tested this value by carrying out an unweighted 
least-square fit to the mean radial velocities of a UMi 
A listed in Table II of Kamper (1996). In this solu- 
tion we solved for T pcr i only, while we adopted all the 
other spectroscopic elements as given by Kamper (1996). 
We obtained T pcr i = 1987.63 ± 0.25, in good agreement 
with Eq. (19). However, the formally most accurate ra- 
dial velocity listed in the last line of Kamper's Table II 
does not fit perfectly (O-C = -0.15 km/s) his final or- 



coincidence: (1) Kamper has adopted this value of T por ; as bit with T por ; according to Eq. (19), but rather indicates 

the value of T por ; = 1987.27. The independent radial- 
velocity data published by Dinshaw et al. (1989) lead us 



a fixed input value from Roemer. Nothing is said about 
this in his paper. (2) Kamper found from a full least-square 
solution by chance the same values for T pcr ; and its mean 
error as quoted for Roemer. Such a mere accident is highly 
improbable. (3) The identical values of T per ; and its mean 
error in the two columns of Kamper's Table III occurcd 
due to a mistake or misprint. However, Kamper has not 
published any erratum in this direction. 

Dr. Karl W. Kamper died in 1998 (Bolton 1998). We 
tried to get clarification on the problem of T pcr i from col- 
leagues of Dr. Kamper, but they were unfortunately un- 
able to help us in this respect. Hence we are inclined to 
accept the possibility (1). However, even then there is an 
additional problem with the mean error of T pcr i. Kamper 
has obviously overlooked that Roemer (1965) gave proba- 
ble errors instead of mean errors. Hence the mean error of 
Tperi according to Roemer should read ± 0.12 in Kamper's 
Table III. 



to T pC ri — 1987.57 with a very small formal error. This is 
in good agreement with Eq. (19). Hence we have finally 
adopted T per i as given by Eq. (19). An error of ± 0.13 years 
introduces errors of ±0?3 in i and of ±2?0 in ft, which 
are small compared to the errors in i and ft due to the 
uncertainties in a ph (AP) sin i and in the observed value of 
A/i. 

The values of a a sini, Ka, P, and e given by Kamper 
(1996) in his Table III under DDO+Lick Data are unfor- 
tunately not consistent. If we accept Ka, P, and e, we find 
oa sin i = 2.934 AU, while Kamper gives in his Table III 
2.90 AU. In the text of his paper, Kamper gives 2.9 AU 
for a a sin i. Has he rounded 2.934 to 2.9 and later inserted 
this rounded value as 2.90 into his Table III ? We prefer 
to trust Ka and e, and hence we use for a a sin i the value 
of 2.934 AU (see Sect. 3.1.4) in our investigation. 
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3.2. The astrometric orbit 

3.2.1. Determination of the inclination i 

In Table 3 we compare the observed values of A/i to t (from 
Sect. 3.1.3 and Table 2) with the predicted values of A/xtot 
(using the procedures described in Sect. 3.1.5 and the el- 
ements P, e, Tpori, u> and a p h(AP) sini given in Table 4) for 
different trial values of the inclination i. The length A/z to t 
of the vector A/i is obviously not a function of the nodal di- 
rection Q. The mean error of the predicted value of A/x to t 
includes the uncertainties in all the orbital elements ex- 
cept in i and £1. 

The best agreement between the observed and pre- 
dicted values of A/z to t occurs for i = 50 ?1 (prograde or- 
bit) and for i = 130 ?2 (retrograde orbit). The uncertain- 
ties in the observed value of A/i to t and in the orbital ele- 
ments (mainly in a p h(AP) sini) lead to an uncertainty in i 
of±4?8. 

Our values for the inclination i of the two orbits 
do not fulfill strictly the expected relation i ro trogradc = 
180° ~ ^prograde- The reason is the following: % is determined 
(Table 3) from two slightly different values A/i for the 
prograde and retrograde orbits (Table 2). The difference 
in the A/i values stems from a slight difference in /z and 
hence in /i m (Table 1), and this difference in /io is caused 
by a difference in £ me anph(AP),H (Eq. (3)). The difference 
in the position £ m0 anph(AP),H of the mean photo-center is 
due to the small, but totally different corrections which 
have to be added to the observed HIPPARCOS position 
£ph(AP),av,H m order to obtain the mean photo-center (see 
Table 6 and Fig. 1). As mentioned already at the end of 
Sect. 3.1.1., we had to iterate our procedure of determin- 
ing i and Q, since the corrections depend on these orbital 
elements. 

The fit between the observed and predicted values of 
A/i is rather pleasing. It is not granted that such a fit 
is always possible. In the case of Polaris, for example, 
the spectroscopic orbit and the HIPPARCOS parallax to- 
gether require a minimum value of A/i to t of 2.00 mas/year, 
which occurs for i = 90°. There is no formal upper limit 
for A/itot for i — > 0. However, the requirement that the 
component P is not visible in the combined spectrum 
of AP gave for a main-sequence companion P a spectral 
type later than A8V (Sect. 2.2), or M P < 1.8 M Q . Com- 
bined with the mass function of the spectroscopic orbit, 
f(M) = (M P smi) 3 /(M A + M P ) 2 = 0.02885M©, and 
with a reasonable estimate of Ma(Ma > 5Mq), this 
gives a lower limit for i of about i > 37°, which corre- 
sponds to A/itot < 7 mas/year. Our observed value of 
A/itot of about 5 mas/year fulfills nicely the range condi- 
tion of 2 mas/year < A/i tot < 7 mas/year. 

3.2.2. Determination of the nodal length ft 

Having fixed the inclination i in Sect. 3.2.1, we now deter- 
mine f2 from a comparison of the observed and predicted 



values of the direction ©a^ of the vector A/t. The differ- 
ence (modulo 360°) between the observed value of ©a m 
and the predicted value of ©a^ for ft = gives just that 
desired value of O for which the observed and predicted 
values of @a,x agree. We find = 276 ?2 for the prograde 
orbit and ft = 167 ?1 for the retrograde orbit. The uncer- 
tainties in the observed value of @a^ arid in the orbital 
elements (now mainly in i and T pell ) lead to an uncertainty 
in Q of ±9?5 or ±9?4. 

The quality of the fit in the components of A/i in a* 
and 5 can be judged from the data given in Table 2. The 
overall agreement is quite good. 

3.2.3. The ambiguity problem of i 

If we know only the vector A/t at one epoch and the spec- 
troscopic orbit of a binary, then there is an ambiguity (i 
or 180° — i) in the inclination i, i.e. in the direction of 
motion in the astrometric orbit. In the prograde (or 'di- 
rect') orbit (i < 90°), the position angle of P relative to 
A increases with time, in the retrograde orbit (i > 90°) 
it decreases. The reason for the ambiguity is the fact that 
A/i itself does not indicate whether the orbit will turn to 
the left-hand side or to the right-hand side (see Fig. 1). 

In principle, the knowledge of the mean position of 
the photo-center predicted by the FK5 for T Cj h = 1991.31 
would resolve the ambiguity. However, the mean errors of 
this predicted position of ± 72 mas in a* and ± 66 mas 
in S are so large with respect to the differences between 
%(T c ,h) and i mcan ph(AP)(Tc,H); which are less than 26 
mas (Table 6), that this method is not useful in our case. 

At present, the best solution of the ambiguity prob- 
lem is provided by the results of the photographic obser- 
vations carried out at the Allegheny Observatory, which 
we discussed already in Sect. 2.2 . While the full astro- 
metric orbit based on the Allegheny data (Kamper 1996) 
is not very trustworthy, the Allegheny data give strong 
preference for a retrograde orbit (in contrast to a direct 
one). This can be seen best in Fig. 3 of Kamper (1996): 
The minimum of the residuals (dashed line) occurs for 
i > 120° (cosi < —0.5), and for this range of i the semi- 
major axis derived from the Allegheny data is quite rea- 
sonable. For our preferred value of i (130?2), we read off 
from Kamper's Fig. 3 a value of a p h(AP) ~ 28 mas with 
an estimated uncertainty of ± 9 mas. This is in very good 
agreement with our result, 28.7 ± 2.8 mas. Even the nodal 
length SI derived by Kamper (175°) is compatible with our 
result (167° ± 9°). Kamper's determination of i (179°) is 
very uncertain and therefore not in contradiction to our 
value (130°). He himself says in the text of the paper that 
'all inclinations between 135° and 180° are equally satis- 
factory' in fitting the Allegheny data. (There is a small 
mistake in Kamper's discussion of this point: He claims in 
the text 'that the minimum scatter is for an inclination 
of almost 90°, which results in a face-on orbit'. The rela- 
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Table 4. Orbital elements of a UMi AP 



Table 5. Physical properties of a UMi A and P 



Quantity 


Prograde orbit Retrograde orbit 
(Preferred solution) 


v T = 7 [km/s] 


- 16.42 ± 0.03 


K A [km/s] 


3.72 ± 0.03 


a a sin i [AU] 


2.934 ± 0.028 


a a sin i [mas] 


22.18 ± 1.42 


a p h(AP)sini [mas] 


21.91 ± 1.42 


dph(AP) [mas] 
aA [mas] 

P [years] 
e 

w n 
x n 


28.56 ± 2.73 28.69 ± 2.75 
28.91 ± 2.74 29.04 ± 2.77 

29.59 ± 0.02 
0.608 ± 0.005 
1987.66 ± 0.13 
303.01 ± 0.75 

50.1 ± 4.8 130.2 ± 4.8 
276.2 ± 9.5 167.1 ± 9.4 



tive clause after 90°, his own Fig. 3 and his Table III all 
indicate that '90°' should be replaced by '180°'.) 

A new astrometric space mission will immediately re- 
solve the ambiguity, since it shall then be clear to which 
side of our A/x vector (Fig. 1) the orbit will have turned 
over. Probably the much higher accuracy of a new space 
mission will allow to determine the direction (and amount) 
of the instantaneous acceleration (i.e. to obtain a 'G so- 
lution' in the HIPPARCOS terminology, if not even a full 
orbital 'O solution'). 

3.2.4. Resulting orbits of a UMi AP 

The resulting orbits of the photo-center of a UMi AP are 
listed in Table 4. As explained in Sect. 3.2.3, the retro- 
grade orbit should be preferred. The semi-major axes of 
the orbits of a UMi A and P itself, relative to the center- 
of-mass of AP, and that of P relative to A, are given in 
Table 5. 

In Fig. 1, the two orbits (prograde and retrograde) 
of the photo-center of AP are illustrated. The zero-point 
of the coordinates Aa* and 5 is the HIPPARCOS posi- 
tion xh(T Cj h) at epoch T c .h = 1991.31. The zero-point 
is then comoving with the center-of-mass (cms) of either 
the prograde orbit or the retrograde one. Therefore the 
orbits stay fixed in these coordinates. Since the proper 
motion /i cms of the cms of the two orbits differs slightly 
(Table 6), the linear motion of the position xn(t) predicted 
from the HIPPARCOS Catalogue, differs slightly for the 
two cases. The indicated motion of ih corresponds to A/x 
(Table 2). Hence by construction, the motion of xh is a 
tangent to the corresponding orbit, except for the slight 
difference between the averaged position and the instan- 
taneous position of the photo-center at T c ,h- The dots on 
the orbits mark the positions in intervals of one year, the 
years 1990, 1995, 2000, etc. being accentuated by a larger 



Quantity 


Units 


Combined 


Component 


Component 






A+r 


A 


r 


mv 


[mag] 


+ 1.982 


+ 1.985 


+ 8.5 ± 0.4 


M v 


[mag] 


- 3.63 ± 0.14 


- 3.62 ± 0.14 


+ 2.9 ± 0.4 


M 


[M @ ] 


7.54 ± 0.6 


6.0 ± 0.5 


1.54 ± 0.25 


Spec, type 






F7-F8 Ib-II 


F0V 


Age r 


[years] 




7 TO 7 




a 


[mas] 


142 ± 21 *) 


29.0 ± 2.8 c) 


113 ± 21 c) 


a 


[AU] 


18.8 ± 2.8 *) 


3.84 ± 0.37 c) 


15.0 ± 2.8 c) 


Mean sep. 


[mas] 


135 *) 







orbit of P relative to A; c) orbit with respect to the cms of AP. 



dot. We indicate also the true major axis, on which peri- 
astron, center-of-mass, mean photo-center, and apastron 
are located. In addition we plot the line of nodes. The 
position of the ascending node is indicated by ft. Fig. 1 
demonstrates clearly that the position predicted by HIP- 
PARCOS is drifting away from the actual position of the 
photo-center of AP. 

3.2.5. Derived physical properties of a UMi P 

In Table 5, we summarize some physical properties of the 
components A and P of a UMi. 

The mass of a UMi A is derived from the mass- 
luminosity relation for Ccpheids given by Becker et al. 
(1977). Since we use the luminosity based on the HIP- 
PARCOS distance (132 pc), our value of 6 Mq is higher 
than that of other authors who have used a smaller dis- 
tance. 

The age of a UMi A, and therefore of the whole system 
of Polaris, can be estimated from the period-age relation 
for Cepheids (Becker et al. 1977, Tammann 1969). Using 
Pq = 5.64 days (see Sect. 2.1), we derive an age r of about 
7 • 10 7 years. 

The spectroscopic orbit provides the mass function 
f(M) = 0.02885 M Q . Adopting the inclination i = 130?2 
of the retrograde orbit and Ma — 6.0 Mq for the 
Cepheid, we obtain Mp = 1.54 Mq for the component P. 
Using this value for Mp, we estimate for a star on the zero- 
age main sequence an absolute magnitude of My = +2.9 
and a spectral type of F0V. The magnitude difference 
Amy.AP between A and P is then about 6™5 . As men- 
tioned in Sect. 2.2, a White Dwarf is ruled out by the 
IUE spectra and the low age of Polaris. Our estimate for 
Mp itself would not violate the Chandrasekhar limit for 
White Dwarfs, if we consider the uncertainty in Mp of 
±0.25 Mq. A neutron-star nature of P is possible, but not 
very likely. In any case, the adopted main-sequence nature 
of P is a rather probable solution which is in agreement 
with all observational constraints. Our derived astromet- 
ric orbit does not depend sensitively on the nature of P, 
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since all the possible solutions indicate a very small value 
of (3, so that the difference between the positions of A and 
of the photo-center of AP (see Sect. 3.1.4) is small in any 
case. 

Using Ma and Mp as derived above, the predicted 
semi-major axis of the orbit of P relative to A is 142 ± 21 
mas. Our Table 7 provides a prediction of the position 
of P relative to A, if the ephemerides for Ax orbph (Ap)(i) 
are multiplied by about —4.95. The separation between 
A and P should be presently about 160 mas, is slightly 
increasing to 186 mas until 2006, and is then decreasing to 
about 38 mas in 2017. Hence the next decade is especially 
favourable for resolving the pair a UMi AP. Of course, 
the large magnitude difference of more than 6 m makes a 
direct observation of a UMi P rather difficult. Since A 
and P seem to have nearly the same colour (as judged 
from the spectral types given in Table 5), the magnitude 
difference should be (unfortunately) rather the same in all 
the photometric bands. Nevertheless we hope that modern 
intcrferometric techniques or the use of other devices may 
be able to resolve the pair a UMi AP during the next 
decade. Our paper provides hopefully a fresh impetus for 
such investigations. 

4. Proper motion and position of Polaris 

Center-of-mass of a UMi AP 

The proper motion /i cms (AP) of the center-of-mass (cms) of 
the closest components A and P of a UMi has already been 
derived in Sect. 3.1.1 for the epoch T C , H = 1991.31. This 
proper motion is then transformed to the other epochs by 
using strict formulae, assuming a linear motion of the cms 
of AP in space and time. The values of /i cms (Ap) for the 
epochs 1991.25 and 2000.0 are given in Table 6. 

In order to derive the position a; cms (AP) of the center- 
of-mass of a UMi AP (Table 6), we first transform the 
HIPPARCOS position £ p h(AP),av.H of the photo-center 
of AP from epoch 1991.25 to T CiH = 1991.31 using 
Mph(AP).av,Hi since T c n corresponds best to the effective 
mean epoch of the HIPPARCOS observations. Then we 
subtract from £ p h(AP),av.H(1991.31) the orbital displace- 
ments Ax olh ph (AP).av(1991-31), where Ax is calculated 
from the derived astrometric orbits (prograde and ret- 
rograde), using the averaging method described by Eq. 
(17). This gives us the position x cms (AP)(1991.31) at the 
epoch T Ci h. Using the proper motion Mcms(AP)(^c,H), we 
transform a; cms (AP) 

from the epoch T C , H = 1991.31 to 
the standard epoch 2000.0. For the convenience of those 
users who like to use the HIPPARCOS standard epoch, 
Th = 1991.25, we give also the position x cms (AP) for this 
epoch Tjj. The values which should be used for predicting 
the position a; C ms(AP)(i) and its mean error are given in 
Table 6 in bold face. The right ascension a is given alter- 
natively in the classical notation (h, m, s) and, as done 
in the HIPPARCOS Catalogue, in degrees and decimals 



of degrees. As discussed in Sect. 3.2.3, we propose to use 
preferentially the retrograde orbit. 

The position x C ms(AP)(£) at an arbitrary epoch t can 
be derived by using the strict formulae for epoch transfor- 
mation, using the epochs 2000.0 or 1991.25 as a starting 
epoch. The mean error £ x , cms ,(AP)(i) of a; cms(A p)(t) should 
be derived from 

£ z,cms(AP)(*) = e lcms(AP)( 1991 - 31 ) 

+ 4c ms( AP)('- 1991.31) 2 . (20) 

This equation assumes that Mcms(AP) and 
^cms(AP) (1991.31) are not correlated. This assump- 
tion is not strictly true. However, for most applications 
it is not neccessary to allow for correlations, because for 
epoch differences At = \t — 1991. 31| larger than a few 
years, the second term in Eq. (20) is fully dominating. 
The correlation between /x a *.cms(AP) and ;i{. cms (AP) is 
ncgligably small (only caused by the tiny correlation 
between /io, Q * and jJ,o,s)- 

All the quantities given in Table 6 refer to the HIP- 
PARCOS/ICRS system and to the equinox J2000 (but to 
various epochs). 

4-2. Orbital corrections for the photo-center of a UMi 
AP 

In order to obtain a prediction for the instantaneous 
position Xph(AP)W °f the photo-center of a UMi AP 
at an epoch t, one has to add the orbital correction 
Aa ; orb. P h(AP)(i) to the position of the center-of-mass 

^cms(AP)(0 : 

^ph(AP)(i) = I C ms(AP)(t) + Ax orb , p h(AP)(£) • (21) 

The ephemerides for the orbit of the photo-center of AP 
are given in Table 7. The orbital elements used in calcu- 
lating the ephemerides are those listed in Table 4. Usually 
it is allowed to neglect the effect that the aS system is 
slightly rotating (8 = +0?00088/ycar), due to the motion 
of Polaris on a great circle. Table 7 lists also the position 
of the intantaneous photo-center at periastron, apastron, 
and at T Ci h. The small difference between the instanta- 
neous position and the averaged position (Sect. 3.1.5) of 
the photo-center at T Cj h shows that the deviations of the 
fitting straight line from the actual orbits remain mostly 
below 1 mas within the interval of Dh = 3.1 years of the 
HIPPARCOS observations, since these deviations reach 
their maximum at the borders of -Dh, namely about twice 
the deviation at T Ci h. The very small deviations from a 
straight line explain also why we were, during the HIP- 
PARCOS data reduction, unable to obtain an orbital (O) 
solution or an acceleration (G) solution for Polaris, al- 
though we tried to do so. 

At the end of Table 6, we give the (constant) off-set be- 
tween the mean photo-center and the center-of-mass. All 
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Prograde orbit 



Retrograde orbit 
(Preferred solution) 



Quantity 


Unit 


Epoch 


in q* 


m.e. 


in 8 


m.e. 


in a. 


m.e. 


in 8 


m 


A'cms(AP) 


[mas/year] 


1991.31 


+ 41.17 


0.50 


- 15.49 


0.39 


+ 40.81 


0.50 


- 15.22 


0. 


Mcms(AP) 


[mas/year] 


1991.25 


+ 41.17 




- 15.49 




+ 40.81 




- 15.22 




A*cms(AP) 


[mas/year] 


2000.00 


+ 41.17 




- 15.50 




+ 40.81 




- 15.23 




^phfAPJ.av.H ( + ) 


[mas] 


1991.25 


0.00 


0.39 


0.00 


0.45 


0.00 


0.39 


0.00 





£ph(AP),av,H ( + ) 


[mas] 


1991.31 


0.00 


0.39 


0.00 


0.45 


0.00 


0.39 


0.00 





A£orb,ph(AP),av 


[mas] 


1991.31 


- 10.37 


2.67 


+ 14.89 


3.16 


+ 16.15 


3.08 


- 8.58 


2 


Al mC anpli(AP) 


[mas] 


1991.31 


+ 15.62 


2.62 


+ 12.40 


3.54 


+ 10.62 


3.58 


+ 17.05 


2 


Icms(AP) ( + ) 


[mas] 


1991.31 


+ 10.37 


2.70 


- 14.89 


3.19 


- 16.15 


3.11 


+ 8.58 


2. 


£meanph(AP) ( + ) 


[mas] 


1991.31 


+ 25.99 


2.59 


- 2.49 


4.35 


- 5.53 


4.27 


+ 25.63 


2 



a 



a 



£ph(AP) ,av,H 
#ph(AP),av,H 

^cms(AP) 
^mean ph(AP) 



1991.25 
1991.31 
1991.31 
1991.31 



47 f 075254 
47? 089026 
47? 142856 
47 f 223939 



+ 89° 15' 50 '.'89698 
50 '.'89628 
50 '.'88139 
50 '.'89379 



47f075254 
47?089026 
47?005192 
47^060320 



+ 89° 15' 50 '.'896 
50 '.'896 
50 '.'904 
50 '.'921 



^'cms(AP) 


1991.25 


02 h 


31 m 47 ? 130034 


+ 89° 15' 50 '.'88232 


02 h 


31 m 46 f 992482 


+ 89° 15' 50 '.'905 


^cms(AP) 


2000.00 


02 h 


31 m 48 f 999906 


+ 89° 15' 50 '.'74676 


02 h 


31 m 48 f 846022 


+ 89° 15' 50 '.'772 


^cms(AP) 


1991.25 




37?94637514 


+ 89?26413398 




37°94580201 


+ 89°264140 


^cms(AP) 


2000.00 




37°95416628 


+ 89?26409632 




37?95352509 


+ 89?264103 



Explanation for (+): To the quantities marked with (+) in the second part of Table 6, one has to add the HIPPARCOS positions £ p h(AP),; 
at the corresponding epochs which are given in the first two lines of the third part of Table 6. 



values are valid for the equinox J2000.0, and for the ori- 
entation of the ad system at epoch 1991.31 (which differs 
from that at epoch 2000.0 by AO = -0?008 only). 

The typical mean error of Ax or b jP h(AP)j due to the un- 
certainties in the orbital elements (mainly in f2), is about 
± 5 mas. It varies, of course, with the orbital phase, ap- 
proximately between ± 2 mas and ± 7 mas. However, a de- 
tailed calculation of this mean error is often unnecessary 
for deriving the mean error of the prediction for x ph ( A p) (t) , 
since the mean error of x p h(AP) (t) is governed by the mean 
error of /i cms (AP) for epoch differences \t— T Cj h| larger than 
about 20 years. 

4-3. Comparison of positions 

In Table 8 we compare positions predicted by our results 
with those predicted by HIPPARCOS and by the FK5. 

At epoch T Cj h = 1991.31, the positions of the photo- 
center of AP predicted by our results (for both types of 
orbits) agree with the HIPPARCOS position by construc- 
tion (except for the slight difference between the instan- 
taneous and averaged position). 

From Fig. 1 we see that the HIPPARCOS predic- 
tions for small epoch differences At = \t — T c ,u\, say for 
At < 4 years, are also in good agreement with our predic- 
tions, since the HIPPARCOS data are essentially a tan- 



gent to our astrometric orbits. In other words, the HIP- 
PARCOS data are a good short-term prediction (relative 
to T C) h) in the terminology of Wielen (1997). For larger 
epoch differences (Table 8), the HIPPARCOS prediction 
for x p h(AP)(*) starts to deviate significantly from our pre- 
dictions. Going to the past, e.g. to t = 1900, the differ- 
ences reach large values of about 300 mas = 0'.'3 in each 
coordinate. Such differences are already larger than the 
measuring errors of some meridian circles at that time, 
especially for Polaris. (The formal mean errors in a* and 
S of the position predicted by the linear HIPPARCOS so- 
lution at the epoch 1900 are 43 mas and 50 mas only.) 
The reason for the failure of a linear prediction based di- 
rectly on the HIPPARCOS Catalogue is the fact that the 
quasi-instantaneously measured HIPPARCOS proper mo- 
tion of Polaris contains an orbital motion Afi tot of about 
5 mas/year as a 'cosmic error'. 

Our data reproduce rather well the FK5 positions at 
the central FK5 epochs. This is to be expected, since we 
have made use of these positions in determining ^ ( an d 
hence /U cms ). For T Ci h = 1991.31, the FK5 prediction devi- 
ates rather strongly from our values. This is in accordance 
with the mean error of ^fks and the large epoch differ- 
ences T Ci h — ?c,FK5- (The mean errors in a* and 8 of the 
FK5 position (reduced to the HIPPARCOS system) are 
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Table 7. Orbital corrections Ax or b. p h(AP) (t) for the 
photo-center of a UMi AP 



Table 8. Comparison between predicted positions for the 
photo-center of a UMi AP 



Prograde orbit 



Quantity 
or epoch t 



AXorb, ph(AP)(i) 

1987.00 
1988.00 
1989.00 
1990.00 
1991.00 
1992.00 
1993.00 
1994.00 
1995.00 
1996.00 
1997.00 
1998.00 
1999.00 
2000.00 
2001.00 
2002.00 
2003.00 
2004.00 
2005.00 
2006.00 
2007.00 
2008.00 
2009.00 
2010.00 
2011.00 
2012.00 
2013.00 
2014.00 
2015.00 
2016.00 
2017.00 



Retrograde orbit 
(Preferred solution) 



Prograde orbit 



in a* 



in a* 



Difference 



1987.66 
2002.46 

1991.31 
1991.31 



Periastron 
Apastron 

instantan. 
averaged 



+ 0.44 

- 9.71 

- 14.16 

- 13.97 

- 11.63 

- 8.33 

- 4.61 

- 0.73 
+ 3.18 
+ 7.02 

+ 10.75 
+ 14.31 
+ 17.69 
+ 20.85 
+ 23.77 
+ 26.43 
+ 28.79 
+ 30.85 
+ 32.56 
+ 33.90 
+ 34.81 
+ 35.26 
+ 35.17 
+ 34.47 
+ 33.05 
+ 30.75 
+ 27.37 
+ 22.60 
+ 15.96 
+ 6.90 

- 4.13 

- 6.71 
+ 27.54 

- 10.67 

- 10.37 



- 7.89 

- 3.30 
+ 3.53 
+ 9.49 

+ 14.10 
+ 17.62 
+ 20.27 
+ 22.23 
+ 23.61 
+ 24.51 
+ 24.99 
+ 25.09 
+ 24.86 
+ 24.32 
+ 23.51 
+ 22.43 
+ 21.11 
+ 19.56 
+ 17.79 
+ 15.81 
+ 13.64 
+ 11.29 
+ 8.77 
+ 6.09 
+ 3.29 
+ 0.41 

- 2.48 

- 5.24 

- 7.53 

- 8.57 

- 6.54 

- 5.33 
+ 21.86 



- 7.97 

- 2.18 
+ 5.19 

+ 11.15 
+ 15.51 
+ 18.65 
+ 20.88 
+ 22.39 
+ 23.33 
+ 23.78 
+ 23.83 
+ 23.51 
+ 22.89 
+ 21.98 
+ 20.82 
+ 19.43 
+ 17.83 
+ 16.03 
+ 14.06 
+ 11.92 
+ 9.64 
+ 7.22 
+ 4.70 
+ 2.10 

- 0.55 

- 3.17 

- 5.68 

- 7.89 

- 9.42 

- 9.40 

- 6.08 

- 4.57 
+ 18.73 



15.30 
14.89 



16.60 
16.15 



^^mcanph(AP) 

(1991.31) 



+ 15.62 + 12.40 



+ 10.62 



Epoch 

t 



m a, 



in S 



Retrograde orbit 
(Preferred solution) 
in q* in S 

[masl 



H1PPARCOS prediction minus this paper (instantaneous position) 



- 0.50, 
10.08 
13.71 
12.81 

- 9.93_ 



1900.00 


- 275 


- 330 


- 298 


- 313 


1991.31 


- 


+ 


+ 


- 


2000.00 


- 5 


+ 23 


+ 24 


- 2 


2010.00 


+ 12 


+ 79 


+ 78 


+ 21 


2020.00 


+ 90 


+ 111 


+ 101 


+ 103 



c,a,FK5 — 



1900.00 


- 1 


- 14 


- 2 


- 18 


1916.08 


+ 4 


- 34 


+ 9 


- 42 


1927.17 


+ 8 


- 47 


+ 16 


- 58 


1991.31 


+ 29 


- 126 


+ 60 


- 154 


2000.00 


+ 31 


- 137 


+ 66 


- 167 


2010.00 


+ 35 


- 149 


+ 73 


- 182 


2020.00 


+ 38 


- 162 


+ 80 


- 197 



1900.00 


- 23 


+ 17 


1991.31 


- 1 


+ 1 


2000.00 


+ 29 


- 25 


2010.00 


+ 66 


- 56 


2020.00 


+ 11 


- 8 



- 6.23^K5 minus this paper (mean photo-center): 

- 2.20 
+ 1.91 T 
+ 5.97^, 
+ 9.91 ' 

+ 13.69 
+ 17.26 
+ 20.60 

+ 23.68 

+ 26.5CPrograde minus retrograde orbit (instantaneous positions) 
+ 29.02 
+ 31.23 
+ 33.10 
+ 34.60 
+ 35.69 

+ 36.36 

+ 36.52 
+ 36.13 

+ 35 fif both orbits. The total differences between the prograde 
+ 33.3¥id retrograde orbit are shown at the end of Table 8 for 
+ 30.73>me epochs. An extremum in these differences occurs in 
+ 27.02* (+ 68 mas) and in 5 (- 59 mas) at about the year 2012. 
+ 21.92 
+ 15.03 

_l_ g ^g4- Space velocity of Polaris 

^'From the derived proper motions /i cms (AP) of the center- 

- 7.3if-mass of a UMi AP (Table 6, retrograde orbit), from the 
+ 30.Q-7adial velocity v, = 7 (Table 4), and from the HIPPAR- 

_ g j^jOS parallax ph (Eq. 14), we derive the three components 

- 8.M1 ^! W °f the space velocity v of Polaris (Table 9). We 
neglect a possible intrinsic K term in the pulsating at- 



at the central epochs (given in Table 8) 37 mas and 34 
mas, and at epoch 1991.31 72 mas and 66 mas.) 

How large arc the differences between the positions 
which we predict if we use either the retrograde orbit or 
the prograde one? At T Cj h = 1991.31, the differences are 
nearly zero by construction. At other epochs, the orbital 
differences can be seen in Fig. 1. To these differences in 
the orbital corrections, we have to add the slight posi- 
tional differences which are due to the differences in /it cms 



+ 17 -9nosphere of the Cepheid a UMi A (Wielen 1974). This 

is probably justified, especially in view of the very small 

amplitude of the radial velocity due to pulsation. 

The velocity component U points towards the galac- 
tic center, V in the direction of galactic rotation, and W 
towards the galactic north pole. The velocity vgo is mea- 
sured relative to the Sun. The velocity vlo refers to the 
local standard of rest. For the solar motion we use v = 
(+9, +12, +7) km/s, proposed by Delhaye (1965). The 
velocity vco is the peculiar velocity of Polaris with re- 
spect to the circular velocity at the position of Polaris (see 
Wielen 1974). For the required Oort constants of galactic 
rotation, we adopt A = +14 (km/s)/kpc and B = —12 
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Table 9. Space velocity of the center-of-mass of a UMi 
AP. For detailed explanations see Sect. 4.4 . 





Prograde orbit 




Retrograde orbit 








(Preferred solution) 


Velocity 


U V W 


V 


U V W v 




[km/s] 




[km/s] 


vso 


- 14.4 - 28.2 - 5.5 


32.1 


- 14.2 - 28.0 - 5.4 31.9 


m.e. 


± 1.2 ± 0.8 ± 1.0 




± 1.2 ± 0.8 ± 1.0 


VLO 


- 5.4 - 16.2 + 1.5 


17.1 


- 5.4 - 16.0 + 1.6 17.0, 


vco 


- 8.0 - 16.0 + 1.5 


18.0 


- 8.0 - 15.8 + 1.6 17.8 



(km/s)/kpc. As mentioned in Sect. 2.3, the velocity of the 
center-of-mass of a UMi AP may differ from that of a 
UMi AP+B by a few tenth of a km/s. 

The peculiar velocity vco of Polaris is reasonable for a 
classical Cepheid. According to Wielen (1974), the veloc- 
ity dispersions (<Ju, ay, crw) f° r nearby classical Ccpheids 
are (8, 7, 5) km/s. Hence only the V component of v co of 
Polaris is slightly larger than expected on average. 
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